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Abstract 



In the framework of quantum groups and additive i?-matrices, the fusion procedure al- 
lows to construct higher-dimensional solutions of the Yang-Baxter equation. These solutions 
lead to integrable one-dimensional spin-chain Hamiltonians. Here fusion is shown to gener- 
alize naturally to non-additive i?-niatriccs, which therefore do not have a quantum group 
symmetry. This method is then applied to the generalized Hubbard models. Although the 
resulting integrable models are not as simple as the starting ones, the general structure is 
that of two spin-(s x s') sl{2) models coupled at the free-fermion point. An important issue is 
the probable lack of regular points which give local Hamiltonians. This problem is related to 
the existence of second order zeroes in the unitarity equation, and arises for the XX models 
of higher spins, the building blocks of the Hubbard models. A possible connection between 
some Lax operators L and i?-matrices is noted. 
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1 Introduction 



The construction and diagonalization of integrable one-dimensional spin-chain Hamiltonians 
within the framework of the Quantum Inverse Scattering Method is well-known [|l|, |2|, 0, 0]. 
A given integrable model and all its conserved quantities are encoded in an i?-matrix which 
satisfies the Yang-Baxter equation. The quantum group approach |5|, ^ provides a systematic 
way for obtaining a large class of solutions based on representations of some underlying Lie 
algebra ||7| or super-algebra |8|. By construction, such solutions possess the additivity property 
which means that the initial double spectral parameter dependence reduces to the difference of 
the two parameters. The most famous example is the spin-^ XXZ chain. 

On the other hand models having non-additive ii-matrices have been known to exist for a 
long time. Examples of such models include Shastry's solution for the Hubbard model |^ p!c| ], 
the Chiral Potts models [0, ^ and more recently the Bariev models []l4| , [l^. The 
Hubbard model both in its bosonic and fermionic guises was also generalized to multi-state 
versions while retaining the same algebraic structure. Initial versions were first introduced in 
1 17 1, studied in |p!8| , p!s| , 20 1, further generalized in [21| and fermionized in |^^. (See also |23| for 
another possible fermionization scheme.) All these non-additive solutions of the Yang-Baxter 
equations are isolated and do not yet fit in a general framework. 

A general method for constructing solutions to the Yang-Baxter equation out of a given 
known one is the fusion method. It works by multiplying the same matrix by itself a certain 
number of times, at different values of the spectral parameter, and finally multiplying by a 
projector. This works much the same way as building higher-dimensional representations from 
tensor products of a smaller one and a final projection on a subspace. For instance an sl{2) 
spin-s solution can be obtained by successive fusions of the spin-i solution |24]. It is in fact 
possible to fuse an arbitrary product of R matrices to obtain solutions of the YBE corresponding 
to most (but not always all ||2^ ) representations of a given Lie algebra . In the framework 

of quantum groups, the direct method for finding i?-matrices with a given Lie algebra symmetry, 
and corresponding to a given representation, consists of solving linear equations Q. This method 
and fusion give the same results. 

In |^8| fusion was shown to work for a class of models which retained only some aspects of an 
sl{m) quantum group structure. Higher-dimensional solutions were obtained by fusion, where 
no quantum group symmetry and therefore no direct method existed. 

In this work I derive fusion equations for non-additive i?-matrices. The results of section ^ 
are quite general and require the starting matrix to satisfy only a minimal number of properties. 
The generalized Hubbard models are shown to satisfy these properties. This allows to construct 
higher-spin Hubbard models which appear as two copies of a multi-flavor spin-(s x s') model 
coupled at the 'free-fermion' point. The coupling does not have the simple structure of the 
starting models. The resulting integrable models have non-additive, unitary, i?-matrices but 
appear to lack the usual regularity property, which would allow to obtain local Hamiltonians. 
The source of this lack of regularity is traced back to the spin-s building blocks. These models, for 
s > 1, are not regular but still allow for local mutually commuting quantities through a limiting 
procedure from a generic q value. It is not clear how to implement this for the higher-spin 
Hubbard models. Possible applications to the Bariev and Chiral Potts models are mentioned 
in the conclusion. A possible connection between some Lax matrices and iZ-matrices is also 
proposed. 
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2 Fusion for non-additive i?-matrices 



Given a non-additive solution of the Yang-Baxter equation, it is possible to obtain new solutions 
provided one has a projector point. Expressions for four fused matrices are found along with 
the equations they satisfy. The results of this section are general and hold without reference to 



any particular model. The notation follows closely that of reference [28|. The word fusion here 
is used in the conventional sense, and not in the sense of |2^. For the additive case, Kulish and 
Sklyanin had already realized that only two properties were needed for fusion to be possible. 
The Yang-Baxter equation has to be satisfied and a projector point must exist (p. 108 of [P). 
Consider a non-additive solution i?(Ai,A2) of the Yang-Baxter equation (YBE) 

-Rl2(Al, A2)i?13(Al, A3)i?23(A2, A3) = R23i^2, A3)i?i3(Ai, A3)i?i2 (Ai , A2) (1) 

Additivity means that for a proper choice of parameterization, and after eventual transforma- 
tions such as a gauge (a special similarity transformation on R) Q or twist (a special similarity 
transformation on R) |^ transformation, one can write -Ri2(Ai,A2) = i?i2(Ai — A2). Most 
known solutions of the YBE are additive. This includes in particular all solutions correspond- 
ing to a quantum group symmetry, Uq{G), where Q is any Lie algebra or super-algebra 0, 
Known non-additive solutions include the class of generalized Hubbard models in their bosonic 
and fermionic forms [ 21 , ^ , the chiral Potts models [^] and the Bariev models [0, |l^, 16 1 . 



To implement fusion it is enough that the solution at hand has a projector point. Thus 
consider any solution R of the Yang-Baxter equation (|l|), which becomes proportional to a 
projector at some special values of the spectral parameters pair (Ai, A2). Specifically, define the 
projector vr^^-' through 

fl(A + p,A)^<7(A)vr« (2) 

and let the orthogonal complementary projector. Both projectors are assumed 

to be independent of A, and p is some fixed value characteristic of the i2-matrix at hand. I have 



then verified that the methods used in |25] can be extended to the non-additive setting. This 
yields the following results. 

Let S be the matrix which diagonalizes both projectors. Define two fused matrices, for 
i = l,2, by 

^<12>3(A, A3) = ^rsVg i?i3(A, A3) i?23(A - p, A3) ^^512 (3) 

The matrices (|^ satisfy a YBE where one space is a tensor product of two spaces: 

-^<12>3(A) A3) -R<i2>4(A, A4) i?34(A3, A4) 

= i234(A3, A4) i?^^^2>4(Ai A4) i?^]^2>3(Ai A3) , i = 1,2 (4) 

We have thus obtained new i?-matrices. If di,i = 1,2,3, are the dimensions of the spaces 1, 2 
and 3, and tr(7r'^^)) = d, then after deletion of the vanishing rows and columns, -R<i2>3(A, A3) is 
a dd^ dimensional matrix for i = 1, and {did2 — d)d^ dimensional for i = 2. 

Note also that there is another possible choice of fused matrices obtained by taking the 
right-hand side of the YBE at the projector point (see (4) of ||2^). There is however no essential 
difference with the foregoing choice. 

One can then fuse two matrices -R<i2>3(A, A3) to obtain the matrix -R<i2><34>(A, Z^) defined 
by: ^ 

-R<i2><34>(A, m) = 5*3"/ 45 -R<i2>4(A, ^ - /?) ^<i2>3(A, /^) vrji 534 , z = 1, 2 (5) 
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These matrices have dimensions for i = 1, and {did2 — d)^ for i = 2. They satisfy two 
Yang-Baxter equations {i = 1,2): 

-^<12><34>(^; /^) ^<12>5('^' ^s) -R<34>5(A*' -^s) 

= -^<34>5(M)^5)-R<12>5('^''^5)^<12><34>(-^'^) (6) 
-^<12><34> ('^i t^) -^<12><56> ('^' -^<34><56> (/^' 

= -^<34><56> (/^' ^<12><56> ^) ^<12><34> (7) 

Assume now that the original ii-matrix is regular and unitary, i.e. 

iil2(/X,;U) =C(/U)7'l2 (8) 

i?i2(A,M)i?2i(^,A) = /(A,^)I (9) 

where V is the permutation operator and R21 = "^12-^12^12- The function /(A, /i) is then 
symmetric in its arguments, and c(;u) is some complex, generically non-vanishing function. The 
fused matrices (|5|) inherit the regularity property (P) provided they are correctly normalized, for 
the corresponding c-function not to vanish. This can be achieved by the following normalization. 
Insert a factor of (/(A + p,fi))~^ in the right-hand side of equation @. As R{fj, + p, fi) is a 
non-trivial projector, unitarity (^) implies that /(/Lt + p, p) vanishes for all values of p. The 
normalization just introduced cancels this zero in the numerator and leaves a regular fused- 
matrix. Incidentally, the symmetry of / implies that f{p, p + p) also vanishes. 
The form ^ can be simplified to a more symmetric one: 

^5l2><34>(A, ^J') = J(A + p p) '^12^ '^34^^12 ^4 ^14(A, P - p) R2a{\ - p, p - p) 

X i?13(A,^)i?23(A-p,/x)^S?7rg 512^34 , i = 1,2 (10) 

where the normalization has been included. 

Let 9i/(Ai, A2) denote the derivative with respect to the z*^ slot (i = 1, 2). Taking the limit 
A — > p for the normalized matrices I find: 

4l2><34>(^,/^) = <p)c{p - p) ^^^1^'^ ~ P\ V,,V2, 5r2VgSi2 53-/4^^34 (11) 
^?12><34>(M>/^) = c(/i)c(/i - p) 7^13'P24 S{:^Tlf^Si2 S^lnf^Su (12) 

The unitarity property is inherited independently from the normalization: 

N7?(i) ( ^^ _ /(A - p, ^)/(A, /i)/(A - p,p- p)/(A, /i - p) 

^<12><34>(A,/ij/t:<34><i2>(/^,Aj - fiX + p,p)fip + p,X) 

xS^^^nf^Su 3^141^34 , i = l,2 (13) 

where -R<34><12>(A, /^) = ^13^24 i?<i2><34> (A, P) ^13^24- 

In the framework of the QISM, the quadratic Hamiltonian density of such integrable hierar- 
chies is the derivative at A = /i of the matrix R{X, p) = VR{\, p). Taking the limit yields: 

^;^^<12><34>lA,MjU., - 25i/(;U + p,/x)^<12><34>^^'^^ 

^ 2d,f{p + p^p) ^l"2'^3-/-S4? 1^ (^32(A, /. - /,)i?13(A, p) 

x^24(A - - p)i?23(A - p,/u)) vr[2-'vr^4 5'i25'34 (14) 
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The first term is proportional to the identity in the fused spaces and may be dropped. The 
results of [28| in the additive case can be recovered by setting /(A, fi) — > /(A — fi), with / now 
an even function. (There is an obvious misprint in formula (17) of |2^: /'(O) and /"(O) should 
be replaced by f'{p) and f"{p)-) 

An implicit assumption, which does not affect the Yang-Baxter equations, was made when 



deriving the regularity equations and (12): /(A + /), /i) was taken to vanish like (A — ^u) 
for A — /i. However this zero can be of second order. This will be case for the Hubbard 
models studied in the next section. A zero of any order does not affect the unitarity equation 
(|l^ ) as numerator and denominator compensate to give a finite non- vanishing result. But it is 
necessary to do a second or higher order expansion to find the appropriate expressions for the 
regularity equations. It does not appear possible to prove in general, and with a minimal set 
of assumptions, that regularity still holds. But the result is finite. An argument in favor of 
regularity is the unitarity equation which says ^12(1^,^) squares to the identity. However this 
turns out not to be enough and specific counter-examples are the multi-flavor spin-s models at 
7 = 7r/2, for s > 1. Note that the above issues of regularity will arise for all fused matrices with 
non-simple zeroes in the unitarity equations of the starting i?-matrix. 

In the case of a second order zero for /(A,//), and provided the fused matrix is regular, the 
Hamiltonian density is given by: 

X^24(A - - /))i?23(A - p,/x)) 7rf2'7r^4 5'i25'34 (15) 

Regularity implies that the first term is proportional to the identity and may be dropped. 

The projector property is the only additional ingredient needed to construct fused matrices. 
For the i?-matrices based on Lie algebras, the degeneration of a generically invertible i?-matrix 
to a projector is expected. But for non-additive matrices one has to verify in every case whether 
a projector point exists. Note also that the unitarity property ( [I^ ) for the fused matrix indicates 
that it may have its own projector point. This in turn implies that fusion may be continued to 
another level, or even indefinitely as happens in the quantum group framework. 

The above fusing scheme will now be applied to the generalized Hubbard models. 



3 Hubbard fusion and non-locality 

We first recall the construction of the multi-state or multi-flavor Hubbard models in their bosonic 



form |21]. The connection between the L and R matrices is clarified. Fusion is implemented. 
The connection between double zeroes in the unitarity equation and the lack of regularity is 
discussed on specific examples. 



3.1 A Hubbard primer 

The following 'free-fermions' or XX models are building blocks of the Hubbard models. Let n, 
ni and n2 be three positive integers such that ni +712 = n, and A, B be two disjoint sets whose 
union is the set of basis states of C", with card(^) = ni and card(i?) = n2- Let be a square 
matrix with a one at row a and column (3 and zeroes otherwise. Define 

p(i) = {^^'^ ^ ^^'^ + ^^'^ ^ ^^'^) (16) 



4 



p(2) ^ J2 J2 E^P' ®eP'^ (17) 

a,a'^A /3,/3'e-B 

Latin indices always belong to A while Greek indices belong to B. The complex twist parameter 
X is arbitrary. The free-fermions i?-matrix 

R{X) = P^^^ + P^^^ cosX + P^^hmX (19) 

satisfies the additive Yang-Baxter equation: 

RuiX - /x) i?i3(A) R23{fi) = i?23(/u) RnW RMX - //) (20) 

The interpretation of the multiple-flavors in terms of sl{2) states was done in p^ ]. 

Coupling two commuting copies of the foregoing models gives the Hubbard models. Where 
made explicit, the two copies are denoted by unprimed and primed quantities. Let us stress 
that the copies need not be of the same type. For instance, the 'left' copy can be (ni,n2) while 
the 'right' copy is (ra'^,n2) with n not necessarily equal to n'. {rii = n[ = \ correspond to the 
original Hubbard model.) The twist parameters may also differ. One then defines a multi-flavor 
version of cr^, the conjugation matrix 

C =Y^E'^<^ -Y, (21) 

and a diagonal coupling matrix 

/oo' {h) = cosh I + sinh Co C'^ = exp Q Co C'^ (22) 

The parameter h is related to the spectral parameter A by 

sinh(2/i) = [/sin(2A) (23) 

where U is the coupling constant. One chooses for h{X) the principal branch which vanishes for 
vanishing A or U . The Lax operator at site i is equal to: 

io^ ( A) = /oo' {h)RQi{X) Ro'i' (A) loo' {h) (24) 
Their commutation relations for different spectral parameters at a given site are given by 

R{\i, A2) L (Ai) L (A2) =L (A2) L (Ai) i?(Ai, A2) (25) 

1 2 
where L (Ai) = L(Ai) (g) I, L (A2) = I ® i(A2), and 



i?(Ai,A2) = Ii2ihi)hi{h2) 



-Rl3(Al — A2)i?24(Al — A2) + ' ^ ^ "^^ 



sin(Ai -I- A2) 

X tanh(/li + /l2)i?13(Al + A2)Cii?24(Al + A2)C2 ] /i2(-/ii)/34(-/i2) (26) 

This matrix is non-additive as it is not possible to reduce is spectral parameter dependence to 
Ai — A2. It satisfies the regularity property 

R{Xi,Xi) = Vi3V24 (27) 
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and the unitarity property: 



i?12(Al,A2)ii2l(A2,Al) = cos2(Ai-A2) 

X (cos^(Ai - A2) - cos^(Ai + A2) tanh^(/ii - /i2)) I (28) 

The matrix (26) satisfies the Yang-Baxter equation (||), with Aj and hi related through (p^). 

The Hubbard models have a Lax matrix L which differs from the intertwiner R. One may 
wonder which matrix should be a candidate for fusing, and what is the role of the RLL relation 
as opposed to the RRR one. The matrix R satisfies the symmetric (|l]) while L satisfies the 
asymmetric (p5|). This already singles out the former as the natural object to fuse. Another 
compelling reason is that L and (25) are just special asymmetrical limits of R and (||). Indeed 
one easily finds that: 



i?1234(Al,0) = ^^^/l2(/il)i2l3(Al)i?24(Al)/l2(/il) = ^(12)(34) (Al) (29) 



Setting A3 = in (|1]) gives (^). One also has 

iil234(0, A2) = /34(-/i2)i?13(-A2)ii24(-A2)/34(-/i2) (30) 

where now the coupling of the two copies is made on the quantum spaces rather than the 
auxiliary spaces. Setting Ai = in (jlD gives an RLL relation. The corresponding quadratic 
Hamiltonian and all the other conserved quantities are however essentially the same as for the 
auxiliary-space coupling case. Now the R matrix can be seen as an L matrix with couplings on 
both auxiliary and quantum spaces. The price of this symmetrization is the linear combination 
in (|2^ ) and the loss of additivity. The quadratic Hamiltonian density obtained from R at the 
arbitrary regular point A = /i is given by: 

+ ^13 A? (cosh^ h - C2C4 sinh^ h) 



7^24i^2r(cosh^ h - C1C3 sinh^ h) 
1 
2 



sinh(2/i) (piJ(C2 + C4) + pSiCi + Cs, 



U 



cosh 2h 



(- sin(2/x) sinh(2/i) (^13^2^ 7^24^2? + Vi3pifV2APS'' 



+ 2 sin 2n sinh^ h {pi^pil + V13P2I ) 

+ {p^^^ + cos(2/i)pg + sin(2//)pg ) {p^^l + cos(2/x)^5g) + sin(2^)pg) ; 

where h = h{fi) is given by (|2^) and p^jl = VjkPfkCj, i = 1,2,3 and j,k = l,---,4. The 
indices are interpreted as follows: 1 site-m-unprimed-copy, 2 site-m-primed-copy, 3 — > 
site-(m -|- l)-unprimed-copy, 4 — > site-(m -|- l)-primed-copy. It is only at = where this 
expression reduces to the familiar generalized (bosonic) form of the Hubbard Hamiltonians. 



The / factors in (26) combine into a similarity transformation. It is in fact a special type of 



gauge transformation, and an equivalent i?-matrix is given by 
r(Ai,A2) = i?i3(Ai - A2)i?24(Ai - A2) 



+ ^^^7^^^tanh(/ii + /i2)i?i3(Ai + A2)Cii?24(Ai + A2)C2 (32) 
sm(Ai -I- A2) 
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The corresponding Lax matrix /(A) is given by: 

'■(^' °) = '(^' = RuWR^WIM^h) (33) 

Similarly one finds 

''(O' ^) = ^34(-2/i)i?i3(-A)i224(-A) (34) 

cosh n 

with a coupling on the quantum spaces. The regularity and unitarity properties are satisfied 
without modifications. The equivalent (for periodic boundary conditions) Hamiltonian is sim- 
pler: 

|^r(A,,)|.., = ^i3A?+^24Pi? + ^ ^ ^ (35) 

X {Pi3 + cos(2//)pg^ + sm{2fi)p'i^ ) {p^^i + cos(2^)pg^ + sin(2pt)p2^ ) 

The matrix r shows that the Hubbard structure lies in the linear combination of two objects 
rather than the factors of I{h). 

3.2 Fusion 

The right-hand side of (^) shows that Ai — A2 = ±7r/2 are possible projector points. Both 
values actually yield projectors with the same dimensionality. This result is peculiar to the 
underlying XX system, but there is otherwise no essential difference between the two projectors. 
For definiteness p = +'7r/2 is considered below. Let 



TT 



(1) _ I ( 6(1) , 6(3) N ( p,m , 6,(3) 



1234 



{x + X ^){x' + (x') 1) 



TTTn {P^3' + Pis' ) (^'24 + ^'2/ ) (36) 



^S'L = 1-434 (37) 

The function g{X) defined in is constant and equal to (x + x~^){x' + {x')~^). The expression 
(|3^) is an decoupled product of a projector for each copy of a free-fermion system. To arrive at 
this result one uses the following relation which is proven by a direct calculation: 

[7rg7rS\/i2(/i)/34(-/i)] =0 , V/iGC (38) 

where vr^^^ = (PiP ~^ PiP^ projector of one copy (x may have a different value for 

each copy). 

The dimensions of these projectors are given by their traces. In particular for the unprimed 
copy, one has tr (vr*^^)) = nin2 and tr (tt^^^) = n\ + nin2- The matrices which diagonalize 
one copy of both projectors are given by: 

S = ^ ^'^'^ (g) E"'"' + ^ S^/^ ^ ^/3'/3' 
a,a' 13,0' 

a f3 

(^"^ ^ - ^ p^"^) (39) 

a f3 

a,a' 13,(3' 



+ 
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The diagonalized projectors read 



a,a' I3,P' 



E- 



a /3 



(41) 
(42) 



To use the fusion formulae for the Hubbard models one doubles every space to unprimed and 
primed copies. For instance -Ri4(A, fJ- — p) in ( [l0| ) is replaced by Rii'44/{X, /i — p) or rii/44/(A, ^ — 
p), obtained from (^6|) or (|3^), respectively. The calculations involved in ^ and (|l^) are 

p|j40|) . Similarly, 



straightforward and can be carried out using the explicit expressions (36,37 



the quadratic Hamiltonian density can eventually be obtained from the right-hand side of (14) 
or (^), or directly once ([lO| ) is calculated. The explicit expressions in terms of S-matrices are 
however unwieldy, complicated, and unenlightening unless used for specific applications such 
as writing down the Hamiltonian in terms of higher-spin sl{2) generators, or for an explicit 
diagonalization. 

It is easy to verify that the relation between the fused matrices based on R and r is a 
non-diagonal similarity transformation: 



R 



<12><34> 



(Al,A2) — 5';^2^'S'i'2"S'3/5'3/4/ Al'(^l)-^22'(~^l)-^33'(^2)-^44'(~^2) 
X Si2Sii2'SuSy4' ?^<12><34>(^1; -^2) 'S'i2^'S'i/2"S'34^5'3,^, 
X /ll'( — /ll)/22'(^l)-^33'(~^2)-^44'(^2) 5'i25'i/2' 5*3453/4/ 



(43) 



To unravel the coupling structure of the models just obtained from fusion, we can look for Lax 
matrices as in (29). The following relations for (19) are easily derived: 



i?i2(A±-)Ci 



Ri2{X±-)Ci 



-CiRi2{X T 
-Ri2{±l- 



-) 

2' 

A)C2 



(44) 
(45) 



(Such relations clearly have fermionic counterparts |^2[.) One can then obtain the fused matrices 
( p^ ) (for (p6D) at some particular points: 



<12><34> 



(A,0) 



1 



S12 Sii\iIlV { — h)l22' {h)Si2Sv2' -R<i2><34> W 



cosh h 



-^<l'2'><3'4'>(-^)'^12^'S'j^'2'-^ll'(^)-^22'( — /i)5'l2 5i/2/ 



(46) 



and 



R 



(i) 

<12><34> 



(A, 



vr , 



^^^^2 ^ 5'l2^5'i/2/Al'(/i)^22'(-M5'l25i/2/ -R<12><34>(A " ■^) 



R 



<1'2'><3'4'>(A - ^)5'l2^5'i/2/-^ll/(-/i)/22'(^)'S'l25l/2/ 



(47) 



These expressions correspond to the decoupled product of two copies of multi-flavor spin-1 
{i = 2), or spin-0 (i = 1) models at their free-fermion point 7 = 7r/2. Indeed, contrary to what 
happens in (|29|), the product of 5/5's on the right is the inverse of that on the right, and as such 
the /-matrices implement an innocuous gauge transformation rather than a coupling. (Another 
explanation of the non-coupling nature of the I's is found in the following paragraph.) This 
negative result can be understood with hindsight. A simple coupling through the /-matrices 
would have been naive because the conjugation operator C has very special properties with 
respect to the i?-matrices and is peculiar to the spin-^ representation. 
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The points = 0, ±7r/2, and similarly A = 0, ±7r/2, are then decoupling points. For a generic 
pair (A, fi) there is no decoupling, and (|3|,|^) applied to the Hubbard models yield the multi- 
flavor spin-(0 X i), spin-(l x i), spin-(0 x 0) and spin-(l x 1) Hubbard models. The structure 



of their corresponding matrices has a simple interpretation. As emphasized in section 3.1, it is 



the special linear combination of RR and RCRC appearing in (p3) which is characteristic of 



the Hubbard models. Looking back at (10) and expanding the product with R replaced by r, 
one finds a special linear combination of sixteen terms, each a product of RR and RCRC. The 
sixteen terms precisely account for all possible combinations one could have found natural to 
consider. Both R and RC satisfy Yang-Baxter equations ||lO| , and are therefore candidates 
for linear combinations. Fusion nicely retains this structure. It also shows the preponderant 
role played by the underlying sl{2) structure of the XX models which can be repeatedly fused 
to reach any spin and therefore any spin-(s x s') Hubbard model. 

3.3 Non-locality at = -1 

An important issue is whether a given i?-matrix is regular. Within the QISM, integrability 
ensures the existence of a large number of commuting quantities. The Yang-Baxter equation 
implies that the transfer matrices for sites, t(A,/u) = Tro [RoN{^,fJ-) • ■ ■ ^01(^^5 m)]) mutually 
commute at arbitrary values of either of the spectral parameters (the other one remaining 
fixed). The matrix r(A,/i) is therefore the generator of mutually commuting quantities. Such 
quantities are generically non-local. For periodic boundary conditions, local, commuting spin- 
chain Hamiltonians can however be defined by taking the derivatives of the logarithm of the 
transfer matrix at a point where the i?-matrix is regular. The transfer matrix at such a point is 
proportional to the unit-shift operator on the chain. Its inverse, in the logarithmic derivatives, 
'cancels' most of the operators in the numerators at the sites where no derivative has been taken. 
In the following cases regularity does not hold. Unitarity by itself is not enough to ensure that 
the transfer matrix is invertible, and each case should be considered separately. 

Before turning to the fused Hubbard models, consider the spin-1 matrix which can be ob- 
tained by fusion from the XX models ([l9|) , and which appears in (|6|,|3) for i = 2. The following 



gauge transformation turns the asymmetric m = 2 XXC i?-matrix of |^8|] into a symmetric ma- 
trix 

R^'){\) = (I (g, A{\)) R{\) (I ® A{-X)) (48) 
where ^(A) = Y.a, E'^^'^'e^^^^ + E'^^'^^e'^^^ and ca - ci = 1. At 7 = 7r/2, 

(A) = sin 7 + P(2) sin(7 + X) + P^^) A (49) 

reduces to ([T9[). The net effect on fusing the symmetric version, for any value of q, is to remove 
(before —1) all factors of y^^ and q^^ from (38) in |28]. The symmetric spin-1 matrix is 



reproduced in the appendix. For 7 = 7r/2 the two simple zeroes of /(A) = sin(7 -|- A) sin(7 — A) 
become a double zero at A = tt/2. Setting 7 = 7r/2 in (|5^) allows to cancel out a factor of 
sin A. The resulting matrix satisfies the Yang-Baxter equation, is unitary but not regular at any 
value of A. However one can still define local Hamiltonians through a limiting procedure from 
generic values of 7. One drops the prefactor in the left-hand side of (|5ll), and calculates the 
local conserved quantities with this renormalized i?-matrix: 

i/p+i = (sin7sin27)P— log(Tro[PoJv(A)---i?oi(A)])|;,=o - P>0 (50) 

These commuting local Hamiltonians are finite and non-trivial as 7 — > 7r/2. (The factor 
(sin 7 sin 27)^ may cancel some contributions but leaves the main ones.) Thus despite the lack 
of a regular point it is possible to define local conserved quantities. 
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The function /(Ai, A2) in (p^) has a double zero, at Ai — A2 = ±7r/2, which is inherited from 



the XX models forming the Hubbard matrix. The special cases (^,47) are non-regular matrices 
and indicate that the fused matrix ( p!o| ) (i = 2) for the Hubbard models, is probably not regular 
for all values of A = /i. (For the spin-(0 x 0) case, i = 1, dimensionality considerations imply 
regularity, just as for the XX models.) A definite proof of non-regularity for generic values of 
the spectral parameters would be welcome. There is no known 'quantum' deformation of the 
generalized Hubbard models to invoke a first order zero and take the limit. Although the absence 
of local quantities would be surprising, it is not clear whether they exist and how they can be 
calculated for all the Hubbard models corresponding to coupled spin-(s x s') XX models, with 
s or s' > 1. 

Finally, consider the point = — 1 (7 = ±vr/2) for all the models considered in [^]. The 
function /(A) in the unitarity equation picks up a double zero at A = 7r/2. So one can expect the 
loss of the regularity property for all the matrices -R<i2><34>(-^) obtained from the i2-matrices 
of the defining representations of the multiplicity Am models. This was seen above explicitly 
for the spin-1 Ai model. Further fusions will propagate this non-regularity to all the higher 
representations, with the combined appearance of higher harmonics of 7, that is of higher roots 
of unity for q. Local conserved quantities should however still be obtainable through the limiting 
method described above. That the fourth roots of unity, and more generally n*^ roots of unity, 
play a specific role for the single fiavor models {rii = 1, i = 1, • • • ,m + 1) is not surprising. 
The representation theory of the quantum algebra Uq{sl{m + 1)) is in one-to-one correspondence 
with the one for the undeformed algebra when q is not a root of unity. This is not the case for 
roots of unity: the representation theory becomes richer and more complicated. The defining 
representations are however undeformed for all values of q. 



4 Conclusion 

The fusion method was shown to generalize naturally to non-additive solutions of the Yang- 
Baxter equation. Expressions for the fused matrices, the regularity and unitarity equations 
and the quadratic Hamiltonians were obtained. The issue of non-simple zeroes was raised and 
connected to a possible lack of regularity of the fused matrices. This raised the issue of existence 
of a set of local commuting quantities. The generalized Hubbard models were then shown to 
allow fusion for all spin-(s x s') representations, and compact expressions were obtained for the 
i?-matrices corresponding to mixed spin-0, ^, 1 multi-flavor representations. Local Hamiltonians 
are believed to exist but a definite proof and calculation method are lacking. 

The fused Hubbard models inherit the symmetries of the two coupled multi-flavor spin-(s x s') 
copies. These symmetries, and the fusion equations between the various transfer matrices can 
be used to diagonalize them through the algebraic Bethe Ansatz. (The lack of regularity should 
not pose a problem.) 



The connection, noted in section 3.1 for the Hubbard models, between the L and R matrices 
and the corresponding RLL and RRR Yang-Baxter equations could serve as a naturalness test 
for the choice of an i?-matrix. This could be particularly relevant in the Bariev model for which 
more than one i?-matrix is known to exist [14, 15, |T^. More generally one can ask the following 
question. Provided the dimensions match and given a Lax operator associated to a non-additive 
i?-matrix, can L be obtained as a limiting case of the original R or some other one with the 
RLL relation satisfied ? Another general test could be the existence of projector points. Fusion 
should also be applicable to the Bariev and Chiral Potts models. 

Finally, a non-additive matrix for the Bariev model was recently obtained by twisting a 



quantum group i?-matrix and taking a singular limit [30|. It is however not clear whether such 
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a method can be made to work for the Hubbard models. 

Acknowledgement: I would like to thank P. Fendley and I. Horvath for a critical reading of 
the manuscript. 

Appendix: The symmetric multi-flavor spin-1 matrix 



The spin-1 matrix discussed in section 3.3 is given by: 



"^"(^^ - ^) X (X)- (51) 

sin(A + 7) ^<12><34>(A)- (51) 

+ sin(A + 7)sin(A + 27) ^ E"^ E"'^ ® E^" E'^" 

a,b,c,d 

+ sin(A + 7)sin(A + 27) ^ E'''^ ® E^'^ ® E'^'' ® E^^ 
+ sin(27) sin(A + 7) ^ ^ S"'' ® E^'^ ® E""^ 
+ sin(27) sin(A + 7) ^ ^ ® E^"" ® E'"' (g) ^"'^ 
+ sin(27) sin(A + 7) ^ ^ ® ® E''" E^" 

a,b,c a 

+ sin(27) sin(A + 7) ^ ^ ^''^ ® S"'' ® E^^ 
+(sin 7 sin(27) + sin A sin(A + 7)) ^ ^ ^^"^ O 
+ sin 7 sin(27) ^ ^ -E"'' ® ^'^^ (g) ^^'^ 
+ sin 7 sin(27) ^ ^ -E''" ® E'"^ S"'' ® 

a,b a, 13 

+2x-^ cos 7 sin(27) sin A ^ ^ E"" ® E^^" ® E'** (g) 

a, 6 o,/3 

+x sin 7 sin A ^ ^ ^'^^ ^ E'^" ® ^'"^ (g) 

a, 6 o,/3 

sin A sin(A + 7) ^ ^ O E^" 
sin A sin(A + 7) ^ ^ ^"'^ (g (g (g 
sin A sin(A + 7) ^ ^ ® ^^'^ (g (g 
sin A sin(A + 7) ^ ^ ^''^ (g S"'^ (g ^'^^ ® E^" 
sin 7 sin A ^ ^ S"" E^^^ S"'* 

a, 6 a,/9 

+2x^ cos 7 sin(27) sin A ^ ^ E"" ® ® 

a, 6 «,/3 
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+x-^ sin(A - 7) sin A ^ ^ El^f^ E"" ® E^^ 

a,b a,0 

sin(A - 7) sin A ^ ^ E""^ ® E^'' O E"" 

This matrix is regular for 7 7^ ^ + /cvr (A; G Z), and unitary for arbitrary 7: 

^Sg<34> (0) = sin2 7 Pi3^24 vrSJvrg) (52) 
4'iS<34>(A)i??liS<i2>(-A) = sm\j + A) sin2(7 - A) vriJ^f (53) 

where ttC*^) is equal to the right-hand side of (|4^). (It is necessary to let x — > —x before using 
dl) in (||,|7D.) 
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